The Propeller and the Frog 



Margaret Pan 1 and Eugene Chiang 1 ' 2 
mpan@astro . berkeley . edu 

ABSTRACT 

"Propellers" in planetary rings are disturbances in ring material excited by 
moonlets that open only partial gaps. We describe a new type of co-orbital 
resonance that can explain the observed non-Keplerian motions of propellers. 
The resonance is between the moonlet underlying the propeller, and co-orbiting 
ring particles downstream of the moonlet where the gap closes. The moonlet 
librates within the gap about an equilibrium point established by co-orbiting 
material and stabilized by the Coriolis force. In the limit of small libration 
amplitude, the libration period scales linearly with the gap azimuthal width and 
inversely as the square root of the co-orbital mass. The new resonance recalls 
but is distinct from conventional horseshoe and tadpole orbits; we call it the 
"frog" resonance, after the relevant term in equine hoof anatomy. For a ring 
surface density and gap geometry appropriate for the propeller Bleriot in Saturn's 
A ring, our theory predicts a libration period of ~4 years, similar to the ~3.7 
year period over which Bleriot 's orbital longitude is observed to vary. These 
librations should be subtracted from the longitude data before any inferences 
about moonlet migration are made. 

Subject headings: planetary rings, planets and satellites 



1. INTRODUCTION 



Satellites embedded within planetary rings open gaps (IGoldreich fc Tremaind Il980l ). 
Ring particles passing by a satellite are gravitationally repelled so that the satellite is sur- 
rounded by an underdensity of ring material. In Saturn's rings, satellites larger than a 
few km open gaps extending a full 2tt radians in azimuth (e.g., Pan in the Encke gap; 
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Showalter et al.lll986f ). For smaller moonlets, physical collisions and gravitational interac- 
tions between ring particles diffuse particles back into the gap downstream of the moonlet. 
Thus par tial gaps are produced whose azimuthal extents depend on ring vi s cosity and moon- 
let m ass ( ISpahn &: Sremcevial2000l ; ISremcevic et al.ll2002t ISeiB et al.ll2005c iLewis fc Stewart 
2009h . 



T he Cassini spacecraft has imaged such partial g aps (ITiscareno et al.ll2006l ; ISremcevic et al 



20071 ; ITiscareno et al.l l2008r ITiscareno et al.l l2010l . hereafter T10). Because of Keplerian 
shear, density perturbations excited at a moonlet 's position drift toward greater longitudes 
inside the moonlet's orbit and toward lesser longitudes outside, forming a pair of features 
dubbed a "propeller" for its S-like shape. According to numerical simulations, a propeller's 
radial width (the radial o ffset between th e azimuthally extended propeller "blades" ) is 5r ~ 4 
Hill radii of the moonlet ( jSeifi et al.ll2005t ILewis fc Stewartll2009l ) . Combining this result with 
measurements of propellers' radial widths from Cassini images, T10 infer moonlet radii of 
~0.1-1 km. 

The appearance of a given propeller varies significantly from image to image, depending 
on the illumination and viewing geometries, the 3D structure of the density perturbations, 
and the optical properties of ring particles. These effects have not yet been disentangled. 
Still, it seems clear that at least for some propellers, the gaps' azimuthal (longitudinal) 
lengths are much larger than their radial widths. Panel (e) of Figure 1 of T10 shows an 
S-shaped bright lobe of radial width 5r ~ 5.5 km and azimuthal length ~100 km — this 
is the propeller Bleriot, the largest and most extensively imaged — embedded within a still 
longer dark gap whose azimuthal half-length (half the distance between gap ends) is ~ 
500 km ~ 300 Hill radii. 

Intriguingly, Bleriot displays non-Keplerian motion (T10). Over several years, Bleriot's 
orbital longitude deviated from a strictly Keplerian solution, showing residuals that varied 
nearly sinusoidally in time with a period of ~3.7 yr and an amplitude (half the distance from 
peak to trough) of ~0.1 deg or ~200 km, an amount less than but of order L^. Data for 
other propellers are much more limited than for Bleriot but suggest similar behavior (T10). 

The longitude vari ations imply semima j or ax is variations of the underlying moonlet. 
Crida et al.l (120101 ) and iRein fc Papaloizoul ( 120101 ) investigated semimajor axi s variation s 
due to stochastic torques exerted by self-gravitating wakes of ring particles (jSalol 119951 ). 
Over timescales of years, stochastic torques m i ght g ive rise to observable drifts in longitude, 
as shown in Figure 6 of IRein fc Papaloizoul (120101 ). But the observed 3.7-year period of 
the longitude variations does not arise naturally from stochastic torques, which are driven 
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by wakes that have lifetimes of order one orbit period ~ 14 hoursE Alternatively, Bleriot's 
moonlet may librate within a resonance established by another larger moon (T10). However, 
no such resonant partner has been identified. Furthermore, many other propellers show 
longitude deviations of similar magnitude to Bleriot's, and invoking a separate partner moon 
for every propeller seems unnatural. 

Here we propose that Bleriot and propellers like it are indeed participating in resonances. 
But the resonances are with nearby ring material, not with other moons. Each propeller's 
moonlet librates within the potential established by co-orbiting ring particles at the ends 
of the long gaps, at distances away. This new type of co-orbi tal resonance, rem iniscent 
of that between the Saturnian satellites Janus and Epimetheus (jYoder et al. 1983 ). shares 
properties w ith both tadpole and hors eshoe orbits in the restricted three-body problem (see, 
for example, iMurray fc Dermottlll999l ). 



2. THE FROG RESONANCE 

We begin with a toy model of the co-orbital material (Figure [Ih): two identical point 
masses ("secondaries") each of mass /i/2< 1 orbiting a central body ("primary") of mass 
1. The secondaries reside on circular, co-planar, Keplerian orbits of radius 1 and angular 
frequency 1. The azimuthal separation of the secondaries is 2</> <C 2ir. The secondaries are 
meant to represent ring material at the ends of the gap opened by a propeller- moonlet. 

The secondaries establish an equilibrium point between them at azimuth 6 = and 
radius r ~ 1 (in cylindrical coordinates in the frame co-rotating with the secondaries). Just 
as for the triangular Lagrange points, this equilibrium point is dynamically stable because 
of the Coriolis force. And as with conventional tadpole and horseshoe orbits, test particle 
motions about our equilibrium point decompose into a fast epicycle of period 2n and a slow 
libration of the guiding center (semimajor axis). We propose that the slow libration is the 
non-Keplerian motion reported by T10. 

We now derive the libration period and the shape of the guiding center orbit in the limit 
of small libration amplitude. In the frame rotating with angular frequency 1, the equations 
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Fig. 1. — Schematic of our dynamical model for Saturnian propellers. Panel (a) depicts our 
toy model: a test-particle moonlet moving between two co-orbiting point-mass secondaries. 
All lengths and angles are evaluated in the frame co-rotating with the secondaries. Panel (b) 
illustrates how the propeller — the S-shaped figure centered on the moonlet — is embedded 
in a partial gap within Saturn's rings. The point masses in panel (a) represent the ends 
of the heavily shaded, horseshoe-shaped, co-orbital ring in panel (b). The zoomed- in panel 
on the right shows sample trajectories of a test particle computed with the co-orbital mass 
distributed uniformly in azimuth outside the gap (see also Figure Eb) using \i = 10~ 4 and 
L^/r = 0.4. Fast epicyclic motion of different amplitudes is visible on top of the azimuthally 
elongated "frog" orbit of the guiding center. 
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of motion for the (assumed massless) moonlet read 



BU 

r-r9 2 -2r9 = — (1) 
or 

1 dTJ 

r9 + 2r9 + 2r = -— (2) 
r o9 

where U is the celestial mechanician's potential, 

U= 1 - + ^ + ^+ 1 -r\ (3) 
r 2r 2 2r 3 2 ' V ; 

and the distances between the moonlet and the secondaries are 

r 2 = a/1 + r 2 - 2r cos(6> + 0) , r 3 = a/1 + r 2 - 2r cos(# - 0) . (4) 

Note that we are ignoring the displacement of the center-of-mass from the primary. As 
shown below, the associated indirect term of the potential, though crucial for the stability 
of conventional tadpoles, is not essential for our problem. Moreover, in the more realistic 
situation where the mass in the secondaries is spread smoothly over the co-orbital region — 
i.e., over nearly the full 2ir rad in azimuth, excepting the gap of width 20 — the barycentric 
displacement is much smaller than in the present two-point-mass problem. We will consider 
this more realistic case at the end of this section. 

We expand U in the limit of small displacements from the equilibrium point, A = 
r — K 1 and 9 < < n: 



U^ + ^ + !t[i-t±-^- + !L ; ). (5) 



3 3 . 2 /jt / A A 2 9 2 

- + -A 2 + - 1 1 

2 2 V 2 20 2 

We insert (JSJ) into ([I])-®, taking d/dt <^ 1 and keeping only leading-order terms to filter 
out the fast epicyclic motion. We also assume that ^> /i 1//3 , i.e., we assume that the gap 
length is much larger than the Hill sphere of the secondary (not to be confused with the Hill 
sphere of the moonlet). The fixed point is then 

A~|- , 9 = (6) 

00 

and the equations of motion about it are 

9 = -3A/2 (7) 
9 + 2A = (2/i/0 3 )#. (8) 

Equation ([7j) states that the test particle moves according to the Kepler shear. Taking the 
time derivative of (JTJ) and inserting the result into OH]), we have 

= -(6/V0 3 )0 (9) 
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which yields harmonic motion of period 



lib 



— = — 7777 for two point-mass secondaries 
y/3 V 1/2 



(10) 



in units where 2n is the local orbital period. Since A is largest when 9 is largest, the aspect 
ratio of the guiding center orbit is 

max A 4 /i 1//2 



max 9 a/6 4> 3/2 



for two point-mass secondaries. (11) 



Numerical integrations of the two point-mass case agree well with equations ffTUj) (Figure Hk) 
and (fTTj) (data not shown). 

A few comments: 



1. As with conventional horseshoe orbits but not conventional tadpoles, stable librations 
here do not require the indirect potential. 

2. As shown in Figure |2^i, the libration period decreases with increasing libration ampli- 
tude; again, this behavior resembles that of horseshoes but not of tadpoles. 

3. Without the Coriolis force term 2 A in equation (JHJ), stable librations would not occur. 
Thus the Coriolis force acts as a restoring force, just as it does for tadpoles. 

4. Like the triangular Lagrange points, our fixed point corresponds to a local minimum 
of U and a local maximum of the physical potential —U. If the energy of a particle in 
resonance is dissipated slowly compared to the libration period, the libration amplitude 
should grow with time, just as for conventional tadpoles. Our numerical simulations 
confirm this. 

5. Given <p ^> l^ 1 ^ 3 , the guiding center orbit is thinner in r than in 9. Conventional tadpole 
orbits are even thinner: (max A/ max^) ta dpoie = \/37t- 

Given these commonalities with tadpoles and horseshoes, we refer to our guiding center 
orbits as "frog orbits. "□ 

We now address the more realistic situation where the secondary masses are spread 
uniformly in azimuth to fill the region 9 = <ft to 9 = 2tt — <p (passing through 71"). In this 
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called a frog. See, e.g., 



http : //www.horsemanmagazine . com/wp- content /uploads/2008/09/horse-hoof-f rog. jpg. 



-7- 



1— i— |— i— i— i— i— i— i— i— i— i— |— i— i— i— i— i— i— i— i— r 

(o) Two point-moss secondaries 



0.3 - 



- + ^=10" 
A fj.= 10~ 
□ /i=10" 



0.1 - 



', = 0.4 
= 0.4 
= 0.4 

■ ■ ■ I 



^=10 , 0=0.04 
/i =10" 9 1 0=0.04 
+ M=10" 11 , = 0.004 
At =10" 12 , = 0.004 

I ■ 



4, 



-i 

□ 



0/4 0/2 30/4 

azimuthal libration amplitude 



100 



kilometers at the orbit of Bleriot 

200 300 400 



500 



3000 - 



o 



2000 - 



c 
o 
o 



1000 



A fj.= 1 
O m= 1 
* /i=10 
+ /i=10 
X ^=10 



A 



0=0.4 
= 0.4 
= 0.4 
0=4.10" 
= 0.04 
■ ■ ■ I 



(b) Filled co — orbital region - <5 



o 



_L 



_L 



o 



0.000 



0.001 0.002 0.003 

azimuthal libration amplitude»(0. 004/0) (radians) 



— 4 ® 

c 

o 
o 



1- 



- i >- 



0.004 



Fig. 2. — Libration period of a test particle as a function of azimuthal libration amplitude, 
computed from direct numerical integrations for a variety of co-orbital masses /i and gap 
azimuthal half-widths 0. (a) Results for the toy model with two point-mass secondaries. The 
y-axis is scaled to \i and according to equation ( flUj) . whose validity is confirmed by the data 
for small libration amplitude, (b) Results for the smoothed-out case where the co-orbital 
mass fi is modeled as a line source of uniform mass density extending in azimuth from to 
2n — (passing through 7r; see Figured]). The y-axis is scaled according to equation (fl2|) ; 
that all data fall on the same curve confirms the scalings of this equation. Note that in both 
the two-point-mass and smoothed-out models, libration periods decrease as the libration 
amplitude increases. This behavior is similar to that of conventional horseshoe orbits. 
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case we might expect only the mass within an angle ~0 of either end of the gap to influence 
the moonlet. In other words, the effective mass of the co-orbital region may be contained 
at azimuths \6\ between and ~20. If so, then the scalings in equations (ITU]) and ( II ip still 
apply with \x replaced by the effective mass ~(0/7r) ■ fi, where now \x is the mass of the entire 
ring co-orbiting with the moonlet. Our numerical simulations (Figures [lb, |2b) confirm the 
resulting scalings of period and aspect ratio with /x and and yield coefficients such that 

Piib * 6.4 ^ (12) 

— ~ 0.7 — — for fj, = mass of entire co-orbital ring (13) 

max <p 

in the limit of small libration amplitude. For parameters inspired by Bleriot and its environment- 
i.e., for = L^/r ~ 500 km/135000 km ~ 0.00 4 and u ~ 27rSrjr/ M Satura ~ 2.5 x 10~ 12 , 



where £ ~ 30 gem" 2 is the disk surface density (IColwell et al.l 120091 ) and Sr ~ 5.5 km — our 
numerical integrations yield Pub ~ 4 yr. This is encouragingly close to the ~3.7 yr period 
measured by T10. Note that these parameters satisfy the assumptions /i C 1, /i 1 / 3 <(^<1, 
and A < 5r <C 1 that we made to derive equations (IT0l)-(lT3l). Also, assuming a moonlet 
density of ~1 g/cc and the median moonlet radius ~0.75 km inferred by T10, the moonlet- 
to-Saturn mass ratio /i m00 n is ~3 x 10~ 15 <C \x as required by our assumption that the moonlet 
mass is negligible. 



3. SUMMARY AND DISCUSSION 

In this work we identified a new kind of orbital resonance between a moonlet embedded 
within an underdense gap of angular extent <p ^ it and disk material co-orbiting with the 
moonlet beyond the ends of the gap. We found formulae for the resonant libration period 
Pub (equation H2J and the radial/azimuthal aspect ratio of the guiding center orbit (equation 
IT3"]) in terms of and the co-orbital disk mass /i. These librations resemble standard tadpole 
orbits because oc l//i 1//2 and because the Coriolis force stabilizes the resonance. But 
whereas tadpoles also require the indirect potential for stability, the librations discussed 
here do not and are more like horseshoes in this respect. We propose that our new "frog 
resonance" explains the non-Keplerian motion seen for Saturnian propeller features such 
as Bleriot. The observed timescale over which Bleriot 's orbital longitude residuals vary 
(~3.7 yr) is well reproduced by our theory. Measurements of non-Keplerian motion and 
their interpretation as resonant frog librations thus offer new constraints on ring surface 
density (through fi) and gap geometry (through 0) independent of radiative transfer models 
of scattered light images. 
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Our simple resonance model does not include interactions with other satellites, nor does 
it account for disk torq ues responsi ble for Type I or Type II migration (for migration in the 



context of gas disk s, see [W ard 



discussed by, e.g., ICrida et al 



19971; analogous effects for particle disks like Saturn's rings are 
Since no propellers are seen at the locations of strong 



mean-motion resonances with the major Saturnian moons, and since Bleriot at least has no 
obvious mean-motion resonant partner, we believe it reasonable to neglect such resonances. 
Secular interactions do not alter the semimajor axis of a propeller- moonlet; moreover, the 
eccentricities involved are too small for secular precession to be relevant. 

Type I migration is driven by imbalanced Lindblad torques from the disks exterior and 
interior to the moonlet 's orbit. The torque from each side is dominated by material at the 
edge of the gap, located at a radial distance ±5r away from the moonlet. If we assume that 
the torque from the outer disk exceeds that from the inner disk by of order Sr/r, where r is 
the moonlet 's orbital radius, then the Type I migration rate is given by 



r Type I ~ — Mdisk ' Me 



r y 

Sr) 



""Kepler ~ ~1 Cm/ yr, 



(14) 



where //disk = Sr 2 /Ms a turn; £ is the unperturbed local surface density; /i m00 n is the mass 
ratio of the moonlet to Saturn; and f Kepler is the Keplerian orbital veloci ty of the moon - 



let. Equation (fT4j) follows from the standard impulse approximation (e.g., iDermottl 119841 ) . 
Because the gap radial width 5r ~ Mm oonr for propellers, the scalings in equation (Ti~4"|) are 
equivalent to those of equation (39) of ICrida et al.l (120101 ): ^Typei oc /idisk jttmoon ^Kepler- Our 
numerical estimate for rT yP ei in equation ( fl4"l) is made for parameters appropriate to Bleriot 
and is much too small to contribute significantly to that propeller's longitude residuals, as 
those residuals imply an rms radial speed of ~100 m/yr. Moreover the Type I drift is of one 
sign, whereas the radial drift of Bleriot implied by the observations switches sign. 

However, since 5r/r ~ 4 x 10~ 5 , the near cancellation of Lindblad torques assumed 
above is delicate and may be overwhelmed by radial surface density gradients. Any long- 
term trend in the longitude residuals, over and above the sinusoidal frog oscillations which we 
h aye analyzed, may there fore help constrain the radial density profile. For example, Figure 4a 
of iTiscareno et al.l ( 120101 ) shows an apparent monotonic drift in Bleriot 's longitude residuals 
of about +0?05 in 3.7 years. If this drift is real and if it grows quadratically with time — and 
these are highly uncertain prospects, as any characterization of a long-term trend depends on 
the method used to derive Bleriot's mean motion (M. Tiscareno, personal communication)— 
then such a drift would correspond to an inward migration speed of ~3 m/yr. The disk just 
outside Bleriot's orbit would have a greater surface density than that inside by of order 1 
part in 100. Alternatively, a l ong-term nonlinear tren d in the longitude c o uld a rise from the 
stochastic torques studied by ICrida et al.l (120101 ) and lRein fc Papaloizoul (120101 ). 
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We have interpreted the non-Keplerian motions of propeller-moonlets as backreactions of 
their perturbed disks on the moonlets. Our model, however, is not self-consistent because we 
have not considered how the motions of the moonlet feed back into shaping the gap. We have 
only computed the motion of a moonlet, idealized as a test particle, embedded in a gap whose 
structure is presumed stationary in the co-rotating frame. We do not have a theory that 
predicts the libration amplitude relative to the gap ends. Energy dissipation (for example, by 
interparticle collisions) tends to increase the libration amplitude, but we do not understand 
how the amplitude might be damped. Observationally, whether the gap itself shows non- 
Keplerian motion, or equivalently whether the moonlet moves relative to the gap ends, is 
difficult to determine (M. Tiscareno, personal communication). The appearance of a gap 
varies from image to image in ways not yet completely understood, and even measurements 
of basic quantities such as a gap's angular size are complicated by azimuthal structure 
inside the gap (i.e., gap ends are not simple step functions). Improvements in radiative 
transfer modeling of images; more astrometric measurements of propellers and their relative 
positions inside their gaps; and a self-consistent theoretical treatment of how moonlets force 
disks and disks force moonlets may yield further insights. 

We thank Matt Tiscareno for sharing his discovery in advance of publication that pro- 
pellers exhibit non-Keplerian motion, and for subsequent discussions. MP appreciates the 
hospitality of the Canadian Institute for Theoretical Astrophysics, where part of the writing 
was done. We thank Glen Stewart for a helpful referee's report, and Joe Burns for encour- 
aging remarks. We thank Aurelien Crida, John Papaloizou, and Hanno Rein for feedback 
on our manuscript. This research was supported in part by the National Science Foundation 
and UC Berkeley's Center for Integrative Planetary Science. 
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